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Introduction. This paper is devoted to the development of a method for estimating the current time-frequency 
parameters of each of a set of simultaneously and independently operating generators in the radio electronic system. A 
general case is considered, in which the deviation of true values of the nominal generator parameters from the assumed 
values of these parameters is determined not only by random, but also by long-term frequency deviation. The work 
objective is to generalize the method for estimating the time-frequency parameters of signals (long-term nominal 
frequency and current frequency deviation from the nominal value) based on the simultaneous measurement of the 
phases of signals generated in the system of independently functioning generators. The research task is to consider a 
system of simultaneously and independently operating generators. Each of them generates harmonic signals, whose 
time-frequency parameters, such as the average frequency, are constant during a certain interval of observation. But 
herewith, these time-frequency parameters are known with insufficient accuracy due to the influence of external factors 
(changes in temperature, pressure, supply voltage, etc.). It is required to obtain estimates of the time-frequency 
parameters of signals (the duration of the measuring interval, values of the long-term frequency and the standard 
deviation) from the results of measurements of the phases of signals formed by the generators at measuring intervals 
belonging to the observation interval, within which the average frequency remains constant. 

Materials and Methods. A system of simultaneously and independently functioning generators is considered. The long- 
term value of signal frequency for each of the generators over the observation interval remains constant, but it is known 
with some margin of error. During the observation interval, several measurements of the signal phase of each of the 
generators are performed. At the same time, the current values of the signal frequency and the duration of the measuring 
interval have random deviations from the long-term values, and follow the normal distribution law with zero 
mathematical expectation and a known variance. The estimation of time-frequency parameters based on the results of 
measuring the signal phases is carried out using a multidimensional likelihood function. The maximum is found on the 
base of solving the redefined system of linear algebraic equations. 

Results. A new mathematical model and a numerical-analytical method for determining the time-frequency parameters 
of signals are developed. They take into account both the long-term constant frequency deviation and short-term 
random deviations. 

Discussion and Conclusions. The results obtained can be used under the development and creation of data-measuring 
and information-telecommunication systems, including geographically distributed systems. The resulting estimates of 
the time-frequency parameters enable to increase the signal frequency stability and, accordingly, to improve the 
accuracy of measurements and the quality of information transfer. 
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Metox_ OeHHBaHHA § JIMTeIbHOCTH BpeMeHHOTO HHTepBasia Cc  HCHOIb30BaHHeM 
OAHOBpeMecHHBIX H3MepeHHit a3bl B CHCTeMe OAHOBPeMeHHO HM He3aBHcCHMO padoTaolnx 
reHepaTopoB 


A. Al. Ta6pnoubsan, A. A. Kocroraoros’, O. A. Cacapban’, JI. B. Uepkecoga’, O. B. Jisopunkos® 

' @YTI «Pocroscknii-Ha-JJoHy Hay4HO-HccejOBaTeECKHH MHCTHTYT paqMocBs3M» (r. Poctos-Ha-Jlouy, Poccuiickas Denepauua) 
> PP BOY BO «Jlonckoii rocyzapcTseHHblii TexXHu4ecKHii yHuBepcnteT» (r. Pocros-Ha-Jouy, Poccuiickas Pexepaina) 
3OAO «Muncxuii HaydHo-nccneqoBaTenbcKuit MpHOopocTponTenbHBIii MACTHTYT» (r. Munck, Pecrry6m«a benapycp) 


Beedenue. CraTba TocBaleHa pa3paOoTKe MeTOJa OL|CHHBAHHA TEKYLIMX YaCTOTHO—BPeMECHHBIX MapaMeTpoOB Kak Oro 
3 COBOKYMHOCTH OJHOBpeMeHHO HM He3aBHCHMO (byHKUMOHHpyIOWMx TeHepaTOpoB B COCTaBe paj{MOIIeKTPpOHHOH 
cuctembl. PaccmotpeH oOmMit ciyyali, IpH KOTOPOM OTKJIOHCHHe HCTHHHbIX 3HAYCHHM HOMMHAJIbHbIX MapaMeTpoB 
reHepaTOpoB OT ipeqMONaraeMbIx 3HAaYeCHHH yKa3aHHbIX MWapaMeTpOB oOMpesesAeTcA He TOKO Cy4alHbIM, HO U 
JOTOBPeMeHHBIM OTKIOHEHHeM uacToTEI. Llenb1o paOoTs aBaeTcA OOOOMWeHHe MeTOa OL|CHHBAaHHA “YacTOTHO— 
BPCMCHHBIX MapaMeTPOB CHrHaJIOB (AOTOBpeMeHHOM HOMHHAJIbHOM YaACTOTbI HW TeKyWerO OTKJIOHCHHA 4YACTOTbI OT 
HOMHHaJIBHOrO 3HaYeHHA) Ha OCHOBe OJHOBpeMeHHOrO U3MepeHuA cha3 curHaoB, topMupyeMBIX B CHCcTeMe 
He3aBHCHMO PYyHKIUMOHMpyIouHux reHepatopos. IlocraHoBka 3aya4uu — paccMaTpuBaeTCA CHCTeMa M3 OJHOBPeMeHHO U 
He3aBHCHMO (PyHKUMOHMpyIOWINx TeHepaTopos, KaxKbI 13 KOTOPbIX POPpMUpyeT rapMOHHMYeCKHe CHIHAJIbI, YACTOTHO— 
BPeMeHHbIe MapaMeTpbI KOTOPBIX, TaKWe KaK CpeJHAA 4ACTOTA, ABJIAIOTCA B TeYCHHE HEKOTOPOTO HMHTepBalia 
HaOJOAeHHUA MOCTOAHHBIMH, HO H3BECTHBIMM C HeOCTaTOYHOM TOUHOCTBIO H3-3a BJIMAHHA BHELIHHX (bakTOpoOB 
(H3MeHeHHe TeMilepaTypbl, WaBleHua, HallpsKeHuA MMTaHHA U T.J.). TpeOyeTca 10 pe3yibTaTaM u3MepeHHi (a3 
curHaloB, POpMUpyeMBIX reHepaTopaMH Ha U3MePHTeJIbHBIX HHTepBalax, IpHHaexalluxX HHTepBaly HaOMIOAeCHHA, B 
lipeyjesaxX KOTOPbIX 3HAaYeHHA CpeHel 4YACTOTbI OCTAaIOTCA MOCTOAHHBIMH, MOJYIHTb OWCHKH YaCTOTHO—BpeMeHHBIX 
TapaMeTpoOB CHTHaIOB — JIMTeJIbHOCTH VW3MepHTeIbHOTO UHTepBasla, 3HAaYeHHM OTOBpeMeHHO! UaCTOTHI 
CpeHeKBapaTHYeCKOLO OTKIIOHCHHA YaCTOTEL. 

Mamemamuueckaa modelb u memodvi pewenua. PaccMaTpHBaeTcCd CHCTeMa OJHOBpeMeHHO H He3aBHCHMO 
(PYHKUMOHMpyIoWux reHepatopos. JlonroBppeMeHHoe 3HaYeHHe YaCTOTHI CHIHalla [Id KaxKOTO H3 reHepaTOpoB Ha 
WMHTepBasle HaOMIOACHHA OCTaeTCA MOCTOAHHBIM, HO H3BECTHBIM C HEKOTOPOM MorpemHoctH1o. B TeyeHve HHTepBasia 
HaOMIOJCHHA MpOBOAMTCA HeCKONbKO H3MepeHHi ca3bI CurHaoB KaxKNoro M3 TreHepaTopos. IIpu 9ToM TekyuHe 
3HAYCHHA YACTOTbI CHrHaIa HU JVIMTeIbHOCTH HW3MepHTesIbHOrO UMHTepBalla MMeIOT CJlyyaliHble OTKIOHCHHA OT 
ONTOBPeMeHHBIX 3HAYeCHHM WU MOXYMHAIOTCA HOPMA@JIBHOMY 3aKOHY paciipeyesieHHA C HyJIeCBbIM MaTeMaTHYeCKHM 
OoKMTaHveM VM W3BeCTHOM BeHYMHOM Aucnepcun. 

OueHnMBaHve 4aCTOTHO-BpeMeHHBIX TMapaMeTpOB HO pe3yIbTaTaM H3MepeHHA (a3 CHTHaIOB MpOBOAHTCA Cc 
MCHOb30BaHHeM MHOTOMepHON dyHKuUHU MpaByonoqoous. B KayecTBe OWeHOK cilyyaMHOrO OTKJIOHCHHA YaCTOTHI 
CHTHaIOB H JJIMTeIbHOCTeH WHTepBaIOB H3MepeHHaA, BbIOUpAaIOTCA 3HAYeHHA, WOCTABIIAIOWMe MaKCHMyM (yHKUMH 
lipaByonoyzoous. HaxoxyqeHnve MaKCHMyMa IIpOBOAMTCA Ha OCHOBe peLeHHA MepeompeeseHHOM CHCTeMbI JIMHEMHBIX 
anreOpawuecKHx ypaBHeHHi. 

Pesyaomamei ucciedoeanua. Pa3padoTaHa HOBadt MaTeMaTH4ecKad MOJeJIb HM YHCCHHO—aHasIMTH4eCKH MeTOA 
onpeyeeHua YaCTOTHO—BpeMeCHHBIX TapaMeTpOB CHIHaJIOB, YUHTHIBAIOWIMe Kak JOJTOBpeMeHHOe TOCTOAHHOe 
OTKJIOHEHHE YaACTOTHI, TAK HW KPaTKOBPeCMECHHEIe OTKIIOHEHHE, HOCALIMe CILyuaiHbIi XapakTep. 

O6cyacdenue u 3axmouenun. TlomydeHuble pe3yIbTaTbI MOTyT ObITb MCHONb30BaHbI TIpH paspaOoTKe H co3qaHHHu 
MHOpMalHOHHO-M3MePpHTeIbHBIX HW HH(POPMAaI[MOHHO-TesIeKOMMYHMKalHOHHBIX CHCTeEM, B TOM 4uciIe 
TeppHTOpHabHO pacnpeyeueHHEIX cuctem. Ilomyuaemble OLeHKH YaCTOTHO-BPeCMeCHHBIX WapaMeTPOB MO3BOJIAIOT 
HOBBICHTbh CTAaOHWJIbHOCTb YaCTOTbI CHTHaIOB H COOTBETCTBEHHO MOBbICHTb TOYHOCTb MIpOBOAMMBbIX H3MepeHHi 
Ka4ecTBO Mepeqadn HHPopMallHH. 


Kioyveeoie cioea: reHepaTOpbl BbICOKOUAaACTOTHBIX KomeOaHuil, BbICOKOUACTOTHbIC PawHOCHHaJibl cTaTHcTH4eCKHit 
MeTO]I cTaOuuIM3allMu yacTOTBHI, cTaOHJIBHOCTh yYacTOTbl, MCTO], HAMMCHbINHX KBapaToB. 


Alia yumupoeanua: Vadpuonpan, J]. 1. Meroy oneHuBaHua JJIMTesIbHOCTH BPCMCHHOLO HHTepBasla C MCIHOJIb30BaHHeM 
OJHOBPeMEHHBIX H3MepeHHit :ba3bl B CHCTeMe OAHOBPeMeHHO MH He3aBHCuMo paOoTaromjHux reHepatopos // Advanced 
Engineering Research. — 2021. — T. 21, Ne 1. — C. 105-110. https://doi.org/10.23947/2687-1653-2021-21-1-105-110 
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Introduction. The constant growth of requirements for audio and video information transmission, for data 
formation and transmission systems, necessitates an increase in stability of the time-frequency parameters of signals 
generating in radio-electronic systems (RES) [1]. The information-telecommunication systems and data-measuring 
ones are typical examples of such systems. One of the components in providing high efficiency of such systems is 
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associated with the formation of signals with highly stable time-frequency parameters. The latter is relevant, in 
particular, for the systems of information transmission with complex signals, radar and radio navigation systems, 
and audio and video information transmission systems as well. 

Currently, the main approach in the creation of RES is a modular construction principle, which determines the 
use of functionally completed blocks and devices connected to each other. The application of this approach leads to 
presence of a large number of high-frequency signal (HF signal) generators in any radio-electronic system. At the same 
time, functional completeness of each of the devices included in the system enables to consider the HF-signal 
generators of these devices as units operating simultaneously and independently on a certain time interval. 

Despite the fact that signal generators in various devices have different parameters, the presence of a large 
number of simultaneously and independently functioning generators in the RES provides determining the current values 
of the time-frequency parameters of these signals through measuring and subsequent processing of the generated signals 
phases! [2-6]. The latter makes it possible to estimate the current values of the time-frequency parameters with higher 
accuracy. In turn, the obtained values of the time-frequency signal parameters are associated with the parameters of the 
generators that are forming them. That allows either to stabilize the generators frequency, or to consider its current 
value during the subsequent signal processing [2-6]. 

Well-known solutions to the problem of estimating the current time-frequency parameters take into account, as 
a rule, only random components of the frequency deviation due to the influence of various factors. In this case, the 
values of such parameter as the nominal (mean) frequency of the generator are considered known. At the same time, in 
many cases, long-term frequency deviations associated with both the influence of external factors and the technology of 
generator production are not analyzed in the papers under examination. 

Thus, the solution to the problem of estimating the time-frequency parameters of signals, such as the nominal 
(mean) frequency and random deviations from the average frequency of generators, is a challenge. 

The paper objective is the generalization of the method of estimating the time-frequency signal parameters (the 
long-term nominal frequency and the current frequency deviation from the nominal value) based on the simultaneous 
measurement of the signal phases forming in the system of independently functioning generators. 

As in the paper [3], let us consider the system of N +1 generators, each of which forms the harmonic signals, 
whose time-frequency parameters, such as average frequency ,, are constant, but known with insufficient accuracy 


due to the influence of external factors (changes in temperature, pressure, supply voltage, etc.) during a certain 
observation interval. 

It is required to obtain estimates of the time-frequency parameters of signals (the duration of the measuring 
interval, frequency values and RMS (root-mean-square, or standard, deviation) of signals) from the results of 
measurements of the phases of the signals formed by the generators at the measuring intervals belonging to the 
observation interval, within which the average frequency values remain constant. 

Mathematical model and solution method. For the system of N+1generators under consideration, let us 
present the signals formed by each of the generators, as 


s,(t) = A, cos(,-t+9,), n=1,...,N+1. (1) 
Using one of the generators (hereinafter the N+1 generator), we will form measuring intervals of nominal 
duration ¢° (m=1,...,M ), such that during all the specified intervals, the average frequency and RMS of the 


frequency of each of the generators can be considered constant. 


We represent the measured values of signals phase © of each of N generators (n=1,...,N ) on the m-th 


nym 


measuring interval (m=1,...,M ) as follows: 


®,,, =(@, +A,,, )(+At,), n=1..N, m=1,..,M, (2) 


nm 


where Aqw,,,— is random deviation of the n-th generator frequency on the m-th measuring interval from the 


average value w,; f\ and At, are, respectively, nominal duration and deviation from the nominal value of the m-th 


measuring interval duration. 
For the mean value frequency, taking into account the modern production technology of HF—generators, we 


can write 


o,=0 +60,, 60,<o”, n=1,...,N+1, (3) 
(0) 


n 


where the values w\” are known, and 6, , with account for the influence of external factors, are unknown. 


Taking into account the expressions (3), the relation (2) takes the form 


* Gabrielyan DD, Prygunov AG, KhutortsevVV, et al. Generator frequency stabilization method. RF Patent no. 2219654, 2003. (In Russ.) 
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®,,, =(@° +80, +Ao,,, (00 + At,), n=1..,N, m=l,...,.M. (4) 


After linearization (the terms 6@,-At, and Aw,,,-At,, are discarded), the expression (4) enables to estimate 


the random deviation of the frequency of the n-th generator on the m-th measurement interval as follows: 


@, OF 00,4 -o AL, 
he = eee On , n=h.uN, m=1,...M, (>) 


nm 0 
t' ) 





where OD =o -(, 


nm 


The frequency deviations Aw, (n=1,...,.N,m=l,...,M ) and deviations of measurement interval duration 


nm 


At, (m =1,...,M_) follow the normal distribution law [7]: 





2 2 
p(do,,, ) a ex (4®,.) , p(At, )= d ex (4+,.) (6) 


p 2 p > 
(0) (0) (0) (0) o (0) 
270), 2(o “oO ) 21 v4 2(o vat y 





n n m 


where o“” are known values of the relative instability of the frequgncy of the n-th generator (n =1,...,N ). 


The given relations describe the mathematical model of the system of generators that function simultaneously 
and independently. 

Solution method. Taking into account the relations (5) and (6), we can write the multidimensional logarithmic 
likelihood function: 


Z 
1(80, At)= 3° n/m ting, + Lee Pea 80, = “At ) 


> 
m=l n=l 2(o° -@ ao” wy 


n m 


(7) 


which includes the vectors d@ and At, whose elements are unknown values, respectively, 60, (1=1,...,N ) 
and At, (m=1,...,M ). 
The estimates 5, and At, are found from the maximum condition (7) and correspond to the solution of the 


system of linear algebraic equations (SLAE). 











AL (a, At) 
=0, n=l... N, 
050, (8) 
aL (Sa, At) 
=0, m=1,...,M. 
OAt,, 
The equations in (8) have the form 
(®,,,- o” —8e, - tO — ow” -At dit t 
nO nn - “=0, n=1,...,N, 
(0-1) 
(9) 


(® _~o” —60, Ta os a” -At Jo!” 
nym nym m no 0, m=1,...,M. 
(5° -@ o fy 


In a matrix form, the SLAE (9) can be presented as: 
A, A. do _ B, (10) 
A,, A, JL At) |B, )’ 


where 4,, is the block with dimensions N-MxN with elements a"), =(o"” 0); A,, is the block with 





n n 


. . . -2 I 
dimensions N-MxM with elements a‘? =(9"") (a i) : 


mmm 


9 T 
is the block with dimensions N-M x N with elements a‘) = (a‘) : 


nm 


A 


2,1 


is the block with dimensions N-MxM with elements a?” =(o% ee oe ie ; 


mmm m 


A,, 
(0) 


nm nm 


B, is the block with dimensions N-M x1 with elements Db‘? = (69-0) : 
(0) ) (0) 
oa) «tl 


n n 


tO) 
nm QO, m 


B, is the block with dimensions N-M x1 with elements 5 = “(69-2 : 
ot t 


m m 


T is the sign of the matrix transposition. 


Taépuanean J. J. u op. Memoo ouenueanua dnumensuocmu 6pemennozo unmepeana Cc UCnos1b306aHUeM OOHOGPeMeHHbIX UZMepeHUt Pa3zbl 





The system of equations (9) contains N of the unknown components d@,and M of the unknown 
elements At, . The total number of measured values of signal phase is N-M . The estimation of their values, taking into 


account the measurement errors, should be carried out by the least-square method (LSM) [8-10]. In this case, the 
condition N-M >N+M must be met, and the number of measurement intervals must satisfy the next condition: the 
number of measurements M must fulfil the condition M > N/(N-1). 


The representation of SLAE using the expressions (9) and (10) defines the numerical-analytical method for the 
problem solving. All elements of the matrix have an analytical representation. At the same time, when passing to the 
system of normal equations, as a rule, used in the LSM, obtaining the analytical expressions is also not difficult. 
However, the inversion of large-dimensional matrix can be performed only using the numerical methods [11-13]. 


The estimates (Sm,) and (Az, ) obtained from the solution to the system of equations, represented by the 


expressions (9), allow us to determine current time-frequency parameters of the generators and signals forming by 
them. 

Research Results. We propose a mathematical model describing a simultaneous and independent functioning 
of the generator system, and a numerical-analytical method for determining the time-frequency parameters of the 
signals with account for both the long-term constant frequency deviation and the short-term deviation of random nature. 
This approach enables: 

— to evaluate, according to the results of measuring the phases of signals formed by simultaneously and 
independently functioning generators, not only random deviations in the frequency of the generated signals, and to 
obtain the estimates of the average long-term frequency of each of the generators; 

— to exclude an instability influence of the time interval duration of measurements on the resulting 
estimates of the signal time-frequency parameters [14-16]. 

Discussion and Conclusions. The results obtained can be used under the development and creation of data- 
measuring and information-telecommunication systems, including the geographically distributed systems. The resulting 
estimates of the time-frequency parameters provide increasing the signal frequency stability and, accordingly, 
improving the accuracy of the measurements and the quality of information transmission. 
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